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FOREWORD 


The  results  and  conclusions  presented  in  this  report  concerning  the  accuracy, 
response,  and  use  of  the  oil-booted  tourmaline  gage  should  be  of  value  to  those 
interested  in  making  measurements  of  underwater  shock-wave  phenomena.  Useful 
tables  and  curves  for  accuracy  prediction  are  included  along  with  examples  of 
actual  measurements.  This  work  was  funded  through  the  Explosives  Development, 
Effects  and  Safety  Block  of  the  Naval  Sea  Systems  Command  (Task  Area  SF-33-354-391) 
as  a  part  of  the  MADAM  program. 


Approved  by: 

\J.  F.  PROCTOR,  Head 
Energetic  Materials  Division 
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CHAPTER  1 
INTRODUCTION 


There  are  many  considerations  to  be  made  in  selecting  a  transducer  or 
piezoelectric  gage  for  underwater  shock  wave  phenomena  measurements.  The 
accuracy  of  the  total  measurement  system  is  the  primary  factor,  of  which  the 
gage  is  only  a  part.  However,  the  total  system  can  be  no  better  than  the  gage 
or  transducer  selected. 

This  report  addresses  the  characteristics  of  tourmaline  and  tourmaline 
gages,  their  response  and  methods  to  determine  their  response,  the  selection  of 
gage  size  and  predicted  accuracy,  and  the  "real"  shock  wave  and  methods  to 
analyze  and  make  corrections. 
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CHAPTER  2 

THE  TOURMALINE  GAGE 


2.1  CHARACTERISTICS  OF  TOURMALINE 

Tourmaline  is  a  hard  natural  crystal  that  is  bulk  or  hydrostatically 
sensitive. 3-3  Its  piezoelectric  constants  are  high  enough  to  be  useful  and 
are  constant  over  a  wide  enough  temperature  range. 

A  survey  of  piezoelectric  activity  of  minerals  by  W.  L.  Bond^  revealed 
that  only  quartz  and  tourmaline  were  practical  for  extensive  piezoelectric 
work.  Some  of  the  characteristics  for  tourmaline  are: 

Piezoelectric  Const  K  =  11  picocoul/psi 

Dielectric  Const  -  7.44 

Density  -  190-200  lb/ft3  (3.0-3. 2  g/cm3) 

Hardness  -  7.3  (Mohs'  Scale) 

Sound  speed  -  19,685  ft/sec  (6000  m/s) 

Piezoelectric  Moduli  :  d31,  d33 

Tourmaline  must  be  cut  perpendicular  to  the  z  or  optic  axis.  Charge  is 
generated  on  the  faces  perpendicular  to  this  axis  by  pressure  applied  to  these 
faces  as  well  as  the  parallel  faces  or  edges:  K  *  d33  +  2d31.  The  properties 
of  quartz  are  similar  to  tourmaline  except  that  quartz  is  not  hydrostatically 
sensitive,  and  its  edges  must  be  constrained  or  protected  from  the  shock  wave. 
Tourmaline  has  the  hydrostatic  sensitivity  advantage  and  does  not  require  any 
mechanical  shielding  to  permit  a  certain  direction  of  strain,  thus  avoiding  the 
accompanying  problems  of  mechanical  resonances,  reflections,  and  complexity. 


l-Cole,  R.  H.,  Underwater  Explosions  (Princeton  University  Press,  1948). 

2Arons,  A.  B.  and  Cole,  R.  H.,  "Design  and  Use  of  Piezoelectric  Gauges  for 
Measurement  of  Large  Transient  Pressures,"  Cambridge  Thermionic  Corporation 
Publication. 


3Arons,  A.  B.  and  Cole,  R.  H.,  "Design  and  Use  of  Tourmaline  Gauges  for 
Piezoelectric  Measurement  of  Underwater  Explosion  Pressures,"  The  Underwater 

Explosives  Research  Laboratory,  Woods  Hole  Oceanographic  Institution,  NDRC 
Report  No.  A-361,  OSRD  Report  No.  6239,  Mar  1946. 

^Frondel,  C.,  "Tourmaline  Pressure  Gauges,"  The  American  Mineralogist, 

Vol.  33,  Jan-Feb  1948,  pp.  1-17. 
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2.2  CHARACTERISTICS  OF  THE  TOURMALINE  GAGE 

The  bare  crystal  tourmaline  gage  inherently  meets  almost  all  the  desired 
characteristics  for  a  gage  except  for  a  low  impedance  to  match  cables.  The 
impedance  is  high;  it  appears  as  a  voltage  source  in  series  with  a  small 
capacitance.  The  gage  must  be  connected  by  low  noise  shielded  cable  to  an 
impedance  matching  amplifier  and  then  to  more  common  types  of  instrumentation, 
not  a  serious  drawback  in  practice. 

2.2.1  Gage  Coatings 

Many  types  of  gage  coatings  have  been  tested  in  the  past,  and  few  have  been 
found  that  did  not  impair  the  response  of  the  gage  in  some  way.  Some  of  the 
more  serious  impairments  were:  adding  charge  to  the  output,  but  not  in  a 
consistent  way;  slowing  the  rise  of  the  gage  or  response  time;  changing  the 
initial  wave  shape;  absorbing  water  or  leaking,  thus,  lowering  the  gage 
resistance,  the  input  time  constant,  and  impairing  the  low  frequency  response; 
instability  of  gage  calibration  constant  short  term  with  use  and  aging;  and  not 
mechanically  rugged  for  field  test  usage. 

In  the  past,  wax  was  found  to  be  the  best  compromise  coating.  Bostik,* 
some  plastics,  rubbers,  and  epoxies  were  tested  by  several  groups  .^“7  Each 
had  some  of  the  serious  drawbacks  mentioned.  The  wax,  epoxies,  and  plastics 
developed  their  own  charge  in  a  pressure  field.  (Some  of  the  newer  epoxies  show 
promise.)  Reported  tests  with  rubber  showed  slowing  of  the  rise;  RTV's  pass 
water  after  curing  for  some  time  (days). 

The  oil-filled  tygon  plastic  booted  tourmaline  gage  has  proven  to  be  a 
stable,  consistent,  reliable  gage  with  none  of  the  serious  impairments  mentioned 

2.2.2  Gage  Construction 

Gages  for  the  Explosion  Dynamics  Branch  (NSWC)  are  fabricated  from  thin 
circular  discs  of  tourmaline  (See  Figure  1).  The  sizes  more  commonly  used  are: 
1/8  inch  (3.18  mm),  1/4  inch  (6.35  mm),  3/8  inch  (9.53  mm),  1/2  inch  (12.70  mm), 
3/4  inch  (19.05  mm),  1  inch  (25.4  mm).  Gages  up  to  2  inches  in  diameter  have 
been  used  on  occasion.  The  gages  may  consist  of  one,  two,  or  four  piles  of 
discs  with  the  appropriate  tabs  wired  in  parallel.  The  gages  are  usually 
mounted  on  a  nylon  feedthrough  which  is  slipped  into  a  tygon  plastic  tube  or 
boot  filled  with  silicone  oil;  the  closed  end  of  the  boot  may  be  a  molded 
hemisphere  or  a  simple  heated  and  pinched  closure  of  the  cylindrical  boot  (see 
Figure  1).  In  fact,  the  latter  has  proven  to  be  superior  in  reducing 
interference  to  the  shock  wave  and  giving  a  cleaner  response. 


♦Mention  of  proprietary  items  constitutes  neither  endorsement  nor  criticism  by 
NSWC. 

^Arons,  Cole,  "Design  and  Use  of  Piezoelectric  Gauges." 

^Arons,  Cole,  "Design  and  Use  of  Tourmaline  Gauges." 

7jempsey,  J.  B.  and  Price,  R.  S.,  Reduction  of  Scatter  In  Underwater  Shock 
Wave  Measurements  Made  with  PiezoTIectric  Gage,  NOLTR  72-12,  Feb  1972. 


FIGURE  1.  NSWC  TOURMALINE  GAGE 
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2.2.3  Characteristics  of  the  Oil-booted  Gage 

The  gage  calibration  constant  KA(pico-coulombs/psi)  is  a  constant;  it  is 
linear  and  stable  from  about  1  psi  to  5000  psi.  This  has  been  confirmed  by 
calibration  tests  at  NSWC  and  those  reported  in  the  literature.®  It  is 
reasonable  to  assume  that  tourmaline  is  linear  over  a  much  wider  range.  This 
seems  to  be  confirmed  by  the  excellent  agreement  in  experimental  results 
obtained  in  controlled  explosion  tests  at  NSWC. 


Many  gages  have  been  used  and  recalibrated  for  a  span  of  15  years  showing 
no  aging  effects,  with  a  standard  deviation  of  less  than  1  percent. 
Pyro-electric  effects  have  been  reported  to  be  negligible  for  underwater  tests 
for  the  signal  durations  encountered,  even  for  large  charges ,9-il 


The  booted  tourmaline  gage  shows  no  signs  of  hysteresis.  The  upoer 

frequency  response  is  dependent  upon  the  gage  size  or  transit  time  o  he  shock 

wave  across  the  gage,  while  the  low  frequency  is  controlled  by  the  i~  it  time 

constant  of  the  system  to  which  the  gage  is  connected.  With  input  i  dances  of 

100  to  300  Megohms,  and  cable  lengths  normally  required  by  the  safe  >  .Lance 
from  the  charges,  the  input  time  constant  is  100  or  more  times  the  d  ion  of 

the  shock,  usually  yielding  a  lower  response  of  less  than  0.1  HZ.  1  ;e  rise 

or  transit  time  is  considered  in  detail  in  Chapter  3. 


Adequate  sensitivity  with  good  signal-to-noise  ratio  and  an  adequately 
small  transit  time  can  generally  be  achieved.  Sometimes  it  is  necessary  to 
place  a  preamp  near  the  gage  to  meet  these  conditions. 


The  gage's  ruggedness  has  been  proven  in  innumerable  tests,  some  tests  to 
pressures  of  50,000  psi.  The  oil-booted  gage  also  survives  field  handling  well 
with  few  precautions.  The  gage  is  simple  to  construct,  and  advances  in 
electronics  technology  have  simplified  the  impedance  matching. 


Minimum  distortion  includes  a  great  deal.  Any  measuring  device  of  finite 
size  will  disturb  the  shock  medium.  Some  initial  diffraction  of  the  shock  front 
occurs  as  it  passes  over  the  gage.  Some  shock  reflection  occurs  as  the  shock 
front  encounters  the  denser  medium  of  the  tourmaline.  This  also  sets  up 
oscillations  or  multiple  reflections  within  the  gage  discs.  The  Bernoulli  flow 
around  the  gage  causes  the  initial  pressure  to  decrease.  However,  this  has  been 
shown  to  be  a  negligible  effect  for  shock  waves  in  water  up  to  30  Kpsi. 12,13 


®Arons,  Cole,  "Design  and  Use  of  Piezoelectric  Gauges." 
^Cole,  Underwater  Explosions. 

l°Arons,  Cole,  "Design  and  Use  of  Piezoelectric  Gauges." 

^Arons,  Cole,  "Design  and  Use  of  Tourmaline  Gauges." 

^^Arons,  Cole,  "Design  and  Use  of  Piezoelectric  Gauges." 

l^Arons,  Cole,  "Design  and  Use  of  Tourmaline  Gauges." 
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Diffraction  effects  and  internal  gage  reflections  are  considered  in  Sections  3.5 
and  3.6.  Gages  are  used  with  the  disc  edge-on  to  minimize  and  smooth 
distortions.  Face  on  gages  disturb  the  flow,  result  in  internal  multiple  modes 
of  vibration,  and  consequently  do  not  actually  shorten  the  rise  time.  The  rise 
is  more  ragged  and  may  actually  be  longer  for  gages  of  four  discs  face  on.  (This 
is  discussed  in  detail  in  Reference  14.)  Experience  has  confirmed  these 
findings . 


l^Cole,  Underwater  Explosions. 
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CHAPTER  3 
GAGE  RESPONSE 


3.1  THE  SHOCK  WAVE 

The  free  field  underwater  shock  wave  is  generally  described  by  the  simple 
exponential  equation: 


P  *  Pm  e“t/®,  t  >  0  (p  -  0,  t  <  0) 

Where  P  is  the  shock  wave  peak  pressure  at  any  time.  Pm  is  the  peak  shock 
pressure;  t  is  time  from  the  arrival  of  the  shock  front;  and  6  is  the 
shock  wave  time  constant  defined  as  the  time  required  for  the  shock  to 
decay  to  Pm/e,  where  e  *  2.7183. 

This  equation  is  a  fair  approximation  out  to  one  0,  but  for  longer  times  the 
real  shock  wave  decay  does  not  continue  to  decrease  as  rapidly  as  described  by 
this  equation.  This  is  not  important  for  the  upper  frequency  response  of  the 
recording  system  or  to  the  gage  size  selection.  It  does  have  to  be  considered 
in  selecting  the  lower  frequency  limit  of  a  recording  system,  setting  the 
system's  input  time  constant  adequately  high,  and  in  determining  the  error 
introduced  by  the  electronics  system  to  measurements  of  energy  and  impulse. 

3.2  GAGE  SELECTION 

The  selection  of  a  gage  for  underwater  pressure  measurements  generally 
becomes  a  compromise  in  choosing  one  that  is  physically  small,  and  thus  has  a 
short  transit  time  relative  to  the  shock  wave  duration  or  time  constant  6,  and 
yet  large  enough  to  have  an  adequate  gage  constant  (KA)  to  provide  a  recordable 
signal.  Other  factors  also  become  important  in  this  selection:  the  recording 
system's  sensitivity,  noise  levels,  and  frequency  response;  gage  cable  length; 
availablity  of  shock  resistant  preamps  near  the  gage;  the  physical  geometry, 
setup  and  rigging,  charge  size,  range,  cable  signal,  desired  recording  time, 
boundary  and  rig  reflections,  i.e.,  the  total  experiment.  However,  there  is  a 
criterion  for  gage  selection  that  can  be  divorced  from  the  above  factors.  This 
criterion  is  based  on  the  acceptable  accuracy  in  the  measurement  made  by  the 
finite  gage  immersed  in  the  underwater  shock  wave  field. 

3.3  SELECTION  CRITERION 

The  finite  gage  size  affects  the  accuracy  of  measurement  of  the  shock  wave 
peak  pressure  to  the  greatest  degree.  The  peak  pressure  affects  the  measurement 
of  8  directly  (Pm/e).  Impulse  is  the  integrated  area  under  the  curve  to  an 
agreed  upon  point  (typically  3  8  or  10  6),  and  energy  is  the  integrated  area 
under  the  p^  curve  to  the  same  point.  Both  impulse  and  energy  are  rather 
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insensitive  to  small  errors  in  the  shock  wave  peak  pressure.  Therefore,  since 
peak  pressure  is  the  most  sensitive  and  is  an  important  parameter  for  explosives 
comparison  and  effects,  it  is  selected  as  the  measurement  criterion  in  gage  size 
selection.  The  peak  pressure  recorded  is  the  apparent  peak  Papp.  The  actual 
or  real  peak  pressure  is  Pm.  The  response  ratio,  Rp,  is  defined  as  the 
ratio  of  the  Papp  to  Pm.  This  rati'  is  dependent  upon  the  duration  6  of 
the  exponential  shock  wave  and  the  transit  time  of  the  gage,  tD: 

Rp  =  f  (6/tD)  =  Papp/Pm  <  1 


3.4  GEOMETRIC  RESPONSE 

The  thin  circular  discs  that  make  up  the  gage  have  a  geometric  transit 
time,  the  time  for  the  shock  wave  to  cross  the  bare  disc  when  facing  edge  on  to 
the  shock  front  disregarding  all  flow,  interference,  or  distortion  effects  (see 
Figure  A-l).  The  derivation  of  the  equation  for  the  geometric  response  is  given 
in  Appendix  A.  A  Fortran  program  was  written  and  the  computed  responses  were 
obtained.  The  computer  results  consisted  of  the  exponential  plane  wave  response 
of  a  circular  disc  with  50  increments  taken  edge-to-edge  across  the  disc  for 
each  run.  A  separate  run  was  done  for  values  of  0/tj)  ranging  from  1/16  to 
2400  in  small  increments.  The  normalized  or  response  ratio,  Rp,  values  are 
summarized  in  Table  1  for  selected  values  of  S/tp;  Figure  2  shows  the 
plotted  responses  for  selected  values  of  d/tg.  The  rise  with  the  6/tp 
value  of  2400  can  be  considered  the  step  response;  the  differences  in  the  actual 
step  response  values  are  so  small  that  they  cannot  be  plotted  to  the  scale  of 
Figure  2. 

The  response  values  in  Table  1  were  converted  to  percentages  and  plotted  in 
two  ranges  of  e/tp  in  Figures  3  and  4.  The  left  ordinate  is  the  percent  of 
peak  or  X  and  the  right  ordinate  is  the  inverse,  or  the  expected  percent 
error  in  recording  the  peak  of  the  shock  wave  or  %  Error  in  Pk.  Figure  3 
covers  the  more  often  encountered  range: 

5  <  e/tD  <  1000 

Figure  4  covers  the  range  from  0.1  to  10.  The  reduced  time  of  the  peak 
amplitude  is  also  plotted  with  the  t/tp  on  the  right  most  ordinate,  using  the 
same  abscissa  of  0/tp, 

For  further  explanation  of  the  table  and  curves  see  Chapter  4. 

3.5  ACTUAL  MEASURED  RESPONSE 

The  geometric  consideration  is  certainly  a  simplification  of  the  real 
gage.  In  fact,  the  gage  coating  or  covering,  the  oil-boot,  must  be  considered 
as  part  of  the  actual  gage.^  No  covering/coating  has  been  found  that  is 
transparent  to  the  shock  wave.  It  would  have  to  have  the  same  acoustic 
impedance  and  properties  as  the  water,  and  yet  afford  protection  and 
insulation.  The  diffraction  of  the  shock  front  around  the  real  gage  will  also 
tend  to  make  it  appear  larger  or  its  transit  time  longer.  The  strain  and 
signals  induced  in  the  tabs  for  connection  to  the  cable  will  also  make  the  gage 


^Arons,  Cole,  "Design  and  Use  of  Piezoelectric  Gauges." 
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TABLE  1.  ACCURACY  OF  GAGE  RESPONSE 
Response  Rp  vs.  Reduced  Shock  Wave  Duration  0/tD 


e/t  R 

- 2 - E _ 

1/16  .0789 

1/8  .1535 

3/16  .2207 

1/4  .2801 

5/16  .3322 

3/8  .3781 

7/16  .4186 

1/2  .4492 

9/16  .4868 

5/8  .5156 

11/16  .5414 

3/4  .5649 

13/16  .5860 

7/8  .6053 

15/16  .6231 

95999  .6290 

1  .6392 

2  .7880 

3  .8503 

4  .8841 

5  .9059 

6  .9208 

7  .9316 

8  .9398 

9  .9462 

10  .9514 

11  .9557 

12  .9593 

13  .9624 

14  .9650 


e/t 

D 

R 

P 

15 

.9673 

16 

.9693 

17 

.9710 

18 

.9726 

19 

.9740 

20 

.9753 

21 

.9764 

22 

.9775 

23 

.9785 

24 

.9793 

25 

.9802 

26 

.9809 

27 

.9816 

28 

.9822 

29 

.9829 

30 

.983414 

31 

.983940 

32 

.984434 

33 

.984899 

34 

.985336 

35 

.985749 

36 

.986139 

37 

.986508 

38 

.986857 

39 

.987189 

40 

.987505 

41 

.987805 

42 

.988091 

43 

.988364 

44 

.988624 

e/t 

D 

R 

P 

45 

.988873 

46 

.989111 

47 

.989340 

48 

.989558 

49 

.989768 

50 

.989970 

51 

.990163 

52 

.990350 

53 

.990529 

54 

.990702 

55 

.990868 

56 

.991028 

57 

.991183 

58 

.991333 

59 

.991477 

60 

.991617 

61 

.991752 

62 

.991883 

63 

.992010 

64 

.992133 

65  .992252 

66  .992367 

67  .992479 

68  .992588 

69  .992694 

70  .992769 

71  .992896 

72  .992993' 

73  .993087 

74  .993179 
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TABLE  1  (Cont.) 

Response  Rp  vs.  Reduced  Shock  Wave  Duration  6/tp 


e/t 

D 

R 

_ £ _ 

e/t 

D 

R 

P 

e/t 

D 

R 

p 

75 

.993268 

100 

.994923 

240 

.997826 

76 

.993355 

102 

.995020 

250 

.997909 

77 

.993440 

104 

.995114 

300 

.998242 

78 

.993523 

106 

.995204 

350 

.998479 

79 

.993603 

108 

.995291 

360 

.998519 

80 

.993682 

110 

.995374 

378 

.998585 

81 

.993758 

112 

.995455 

400 

.998658 

82 

.993833 

114 

.995533 

450 

.998796 

83 

.993906 

116 

.995608 

480 

.998866 

84 

.993977 

118 

.995681 

500 

.998907 

85 

.994046 

120 

.995751 

600 

.999074 

86 

.994114 

126 

.995949 

700 

.999193 

87 

.994181 

132 

.996128 

720 

.999213 

88 

.994245 

138 

.996292 

800 

.999282 

89 

.994309 

148 

.996442 

900 

.999351 

90 

.994371 

150 

.996580 

960 

.999386 

91 

.994431 

156 

.996708 

1000 

.999407 

92 

.994491 

162 

.996826 

1200 

.999490 

93 

.994549 

168 

.996936 

1440 

.999559 

94 

.994606 

180 

.997134 

1500 

.999573 

95 

.994661 

192 

.997307 

1680 

.999609 

96 

.994716 

200 

.997411 

1920 

.999646 

97 

.994769 

204 

.997460 

2000 

.999657 

98 

.994821 

216 

.997595 

2160 

.999675 

99 

.994872 

228 

.997717 

2400 

.999698 
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look  longer  and  lengthen  the  transit  time.  These  strains  may  also  cause 
transverse  modes  of  vibration  or  oscillation  (more  detail  in  Section  2.6). 

If  the  gage  coating  has  little  deleterious  affect,  then  the  shape  of  the 
rise  of  the  real  gage  should  be  controlled,  or  effectively  be  the  rise  of  the 
bare  crystals  themselves.  Consider  the  real  gage  a  black  box  with  a  transfer 
function  equivalent  to  the  bare  gage.  If  the  actual  signal  recorded  is  the  same 
shape  as  the  signal  predicted  by  computation  except  for  the  duration,  then  the 
coating  has  done  nothing  but  make  the  real  gage  look  like  a  larger  diameter  bare 
gage. 


High  frequency  response  recordings  were  made  with  a  16-Channel  Oscilloscope 
Recording  System  using  fast  sweeps  for  gages  of  size  1/4  inch.  The  shock  peaks 
on  the  film  records  for  the  1/4-inch  gages  are  about  4  cm  in  height.  These  were 
digitized  and  plotted  to  an  expanded  scale  five  times  larger  (20  cm).  The 
1/2-inch  gage  records  were  recorded  on  a  Nicolet  digital  scope  and  plotted  to 
the  same  scale.  Some  records  have  been  obtained  for  1/8-inch  gages  from  large 
shaped  charges  and  tiny  cylinders.  The  few  records  appear  to  have  the  same 
shape  as  those  for  the  1/4-inch  and  1/2-inch  gages,  as  they  should.  However, 
these  have  not  been  plotted  to  the  expanded  scale. 

The  results  of  these  plots  will  be  discussed  in  detail  in  the  next 
section.  It  is  sufficient  to  say  that  the  plotted  real  gage  rise  shapes  were 
excellent  matches  to  the  computed  geometric.  It  should  also  be  noted  that  the 
statistical  approach  used  by  Dempsey  and  Price  in  Reference  16  showed  that  the 
gage  constant  of  the  oil-booted  gage  was  the  same  as  for  the  bare  gage,  and  that 
the  standard  deviation  of  a  set  of  measurements  was  reduced  with  use  of  the 
oil-booted  gage.  Separate  pressure  pot  calibrations  of  the  calibration  constant 
KA  for  the  booted  gages  and  then  their  bare  crystal  elements  were  identical. 

The  rise  responses  of  an  oil-booted  gage  and  a  wax-coated  gage  subjected  to  an 
underwater  shock  are  compared  in  Figure  5.  Notice  that  the  wax  not  only  slows 
the  rise  time  but  changes  the  wave-shape  by  rounding  and  smoothing. 

3.6  RESPONSE  COMPARISONS 

The  measured  transit  time  of  the  oil-booted  gage  is  about  half  again  as 
long  as  the  computed  geometric  transit  time  of  the  bare  gage  element  alone, 
i.e.,  the  time  for  the  shock  wave  to  traverse  the  diameter  of  the  bare  gage 
element:  tp  (measured)  =  1.5  tp  (geometric).  Therefore,  in  order  to 
compare  recorded  gage  response  with  the  computed  geometric  response,  the 
recorded  responses  were  reduced  by  their  respective  transit  times,  i.e.,  an 
actual  gage  rise  was  reduced  by  its  measured  transit  time. 

The  measured  0  was  divided  by  the  measured  tp  to  get  the  values  of  this 
ratio.  The  computer  runs  of  geometric  or  computer  responses  were  identified  by 
this  same  ratio. 

The  comparisons  of  reduced  actual  response  and  geometric  response  were 
plotted  for  matching  values  of  0/tp.  The  vertical  scale  for  pressure 
response  ratio  (defined  in  Section  3.3)  was  also  normalized,  thus,  the  ordinate 

Kpnorm* 


l^Dempsey,  Price,  Reduction  and  Scatter  of  Underwater  Shock. 
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Many  fast  rise  recordings  have  been  made  of  various  size  explosions  and 
compositions.  Typical  charge  weights  were:  1  oz,  1  lb,  5  lb,  10  lb,  15  lb,  50 
lb,  and  1300  lb.  Compositions  were:  C4,  Pentolite,  TNT,  HBX-1,  PBXN103,  and 
Minol.  The  range  of  d/tp  for  those  was  5-90.  A  limited  number  of  fast  rise 
recordings  are  shown  in  Figures  5  through  7  for  1/4-inch  gages.  The  expanded 
and  reduced  plots  for  these  and  a  1/2-inch  gage  digitally  recorded  are  shown  in 
Figures  8  through  11.  Note  how  closely  the  actual  rises  for  the  real  gages 
match  their  appropriate  geometric  curves. 

This  comparison  of  the  actual  rise  shape  to  that  predicted  for  the  bare 
crystals  by  the  geometric  response  proves  to  be  an  excellent  technique  for 
evaluating  gage  coatings  and  various  other  configurations.  This  technique  along 
with  the  comparison  of  calibration  constants  (KA)  for  the  bare  crystals  and  for 
the  covered  gage  defines  the  characteristics  of  the  completed  gage. 

The  small  perturbations  or  oscillations  on  the  plotted  rises  of  Figures  8 
through  10,  are  seen  in  the  actual  oscilloscope  records.  Figure  5  through  7. 
These  are  the  internal  reflections  of  the  shock  front  within  the  tourmaline 
crystals.  The  shock  speed  in  tourmaline  is  about  four  times  faster  than  in 
water.  Note  that  there  are  inflection  points  at  about  the  1/4  points.  The 
reflections  in  the  crystal  would  be  expected  to  be  damped  oscillations  with  the 
same  S-shaped  rise  and  fall. 

Note  also  the  damped  oscillations  after  the  peak.  These  are  at  lower 
frequencies  that  correspond  well  with  the  gage  size  and  the  shock  speed  in  the 
water.  These  are  the  oscillations  mentioned  in  Section  3.5.  They  are  probably 
mechanically  induced  stresses  or  reflections  off  gage  tabs  or  connections. 
Further  investigations  with  modified  gage  constructions  would  shed  more  light  on 
these  conjectures.  Whatever  the  source,  the  oscillations  are  part  of  the  real 
response  of  the  gage;  the  frequency  of  oscillations  is  directly  related  to  the 
crystal  size  and  the  shock  speed  in  water. 
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5  /jsec/cm 


FIGURE  5.  GAGE  RISE  COMPARISON  0009 
OIL-BOOTED  VS.  WAX 


5  fjsec/cm 

FIGURE  6.  GAGE  RISE  SHOT  0012 
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FIGURE  11.  ACTUAL  VS.  GEOMETRIC  RISE 
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CHAPTER  4 

USE  OF  TABLES  AND  CURVES 


4.1  GEOMETRIC  GAGE  RISE  RESPONSES 

Figure  2  shows  the  shapes  of  the  gage  rise  responses  to  one  transit  time; 
it  shows  the  reduced  peak  values  and  the  point  across  the  gage  that  this  maximum 
occurs.  The  gage  rise  responses  shown  are  for  selected  ratios  of  shock  wave 
duration  or  time  constant  0  reduced  by  the  transit  time  of  the  gage.  The 
ordinate  is  actually  the  pressure  response  ratio,  Rp,  the  ratio  of  the 
recorded  or  apparant  peak  (PaDD)  to  the  actual  or  expected  pressure  maximum 

(Pm)- 


Since  the  curves  are  reduced,  they  are  applicable  to  any  tourmaline  gage  of 
similar  construction,  edge  on  to  any  given  exponential  (plane)  shock  wave. 
Chapter  3  showed  that  the  shapes  of  the  rise  for  the  oil-booted  gage  were 
excellent  matches  to  the  computed  geometric  responses.  The  transit  time  of  the 
real  gage  was  about  1.5  times  longer  than  the  computed  transit  time  for  the  bare 
crystal  alone. 

An  example  will  help  clarify  the  procedure  for  use.  Suppose  calculated 
similitude  predictions  indicate  that  0  will  be  63  usee.  A  1/4-inch  gage 
with  a  shock  speed  of  about  5000  ft/sec  would  give  the  geometric  transit  time  of 

4.2  usee.  The  real  gage  will  have  a  transit  time  =  1.5  (4.2)  or  6.3  ysec. 

The  0/tD  would  be  63/6.3  or  about  10.  From  Figure  2,  the  pressure  response 
is  0.95  and  peaks  out  right  at  one  transit  time  or  the  far  edge  of  the  gage. 

This  reveals  that  the  real  gage  should  be  capable  of  recording  95  percent  of  the 
actual  peak  pressure. 

Qualitatively,  Figure  2  reveals  the  detrimental  effect  on  the  rise  shape 
and  peak  values  when  the  gage  becomes  large  relative  to  the  duration  of  the 
shock  wave. 

4.2  ACCURACY  PREDICTIONS  OF  PEAK  MEASUREMENTS 

Figures  3  and  4  show  the  percent  accuracies  and  errors  to  be  expected,  due 
to  the  gage,  in  recording  the  peak  pressures  for  two  ranges  of  shock  time 
constant  (0)  reduced  by  the  transit  time  of  the  gage.  Figure  3  covers  the 
more  used  range  of  5  to  1000,  while  Figure  4  covers  the  range  of  0.1  to  10. 

These  curves  give  the  expected  accuracy  or  errors  in  the  peaks,  but  do  not  give 
any  indication  of  the  shapes  of  the  rises  or  the  times  at  which  they  occur.  The 
curves  do  give  a  quick  way  to  select  a  gage  size  in  the  planning  of  an 
experiment. 
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Take  the  example  in  Section  4.1.  A  shock  duration  of  63  usee  and  the 
real  1/4-inch  gage  response  gives  a  95.15  percent  accurate  peak  measurement. 

But  a  1/2— inch  gage  would  have  a  8/tj)  of  5  and  yield  a  90.45  percent 
accurate  measurement.  This  may  not  be  tolerable.  Now  a  1/8-inch  gage  would 
have  half  the  transit  time  of  the  1/4-inch  gage,  a  6/tD  of  20  would  give 
97.53  percent  accuracy  or  a  2.47  percent  error  in  the  peak  measurement.  A  gage 
selection  based  on  acceptable  accuracy  can  now  be  made;  and  then  one  is  ready  to 
contemplate  the  myriad  of  other  experimental  details,  compromises,  and 
modifications. 

4.3  TIME  OF  PEAK  AMPLITUDE 

The  second  curve  on  Figure  4  is  a  plot  of  the  reduced  time  that  the  peak 
amplitudes  occur  for  reduced  durations  of  O/tp  from  0.1  to  5.  This  curve 
shows  that  the  peak  occurs  at  the  far  edge  of  the  gage  for  values  of  6/tj) 
greater  than  4.  The  ordinate  t/tp  is  at  the  far  right.  These  two  curves 
together  may  prove  helpful  for  unusual  experiments  at  this  extreme,  i.e.,  when 
very  small  charges  must  be  used. 

4.4  RULE  OF  THUMB  ERROR  PREDICTION 

A  handy  rule  of  thumb  to  estimate  the  error  (and  consequently  the  accuracy) 
for  6/ tjj  ratios  larger  than  5: 

%  Pk  Error  *  0.5  _i _ x  100%;  e/tD  >  5 

8/tD 

%  Pk  *  1  -  X  Pk  Error 

In  words,  the  peak  error  prediction  is  one  half  the  inverse  of  the  ratio  of 
shock  duration  0  to  transit  time.  Only  slight  errors  occur  when  this  rule  of 
thumb  is  used.  For  example,  at  «/tD  =  5,  this  rule  predicts  a  10  percent 
error  or  90  percent  of  the  peak.  The  actual  peak  from  Table  1  is  90.59 
percent.  So  the  rule  under  estimates  the  peak  by  0.65  percent.  At  e/tD  = 

10,  the  estimate  is  off  0.15  percent. 

Another  way  to  phrase  this  is  that  the  peak  error  is  the  inverse  of  the 
ratio  0  to  1/2  the  transit  time: 

%  Pk  Error  =  .  1  x  100% 

0/(0. 5  tD) 

The  reason  this  works  is  that  for  0/tj)  greater  than  5,  the  peak  of  the 
rise  occurs  at  one  transit  time  and  falls  upon  the  linear  portion  of  the 
exponential  decay.  An  exponential  decay  is  linear  to  within  0.5  percent  for 
1/10  of  its  time  constant  and  to  within  0.1  percent  for  1/20  of  its  time 
constant.  For  0/tjj  of  5,  1/2  tp  is  1/10  0.  One  half  the  transit  time 
is  the  zero  reference  or  the  center  of  the  gage.  This  procedure  is  analogous  to 
the  lo8 8  in  peak  amplitude  during  the  blank  in  digital  sampling. 
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CHAPTER  5 

ANALYSIS  OF  THE  REAL  SHOCK  WAVE 


5.1  THE  REAL  SHOCK  WAVE 

The  real  shock  wave  is  generally  described  by  the  simple  exponential 
expression  of  Section  3.1.  This  is  certainly  a  convenient  expression  for 
analytical  use,  and  it  is  easily  visualized.  It  is  also  a  fair  approximation 
for  centrally  initiated  spherical  charges. 

The  real  waveform  deviates  from  an  exponential  in  three  ways.  The  two 
commonly  noted  departures  are  at  later  times  near  the  time  constant  6.^ 

The  first  is  that  the  waveform  begins  to  decay  less  rapidly  than  an  exponential 
at  about  6.  The  second  is  the  "hump"  which  is  attributed  to  the  arrival  of 
the  internally  reflected  rarefaction  wave.  The  rarefaction  wave  is  generated  as 
the  shock  front  initially  arrives  at  the  surface  of  the  charge.  (Note  the  small 
hump  at  about  6  in  Figure  7b.) 

The  third  departure  is  at  the  peak  and  initial  decay.  Frequently,  there  is 
a  bunching  and  rounding  or  bending  out,  as  is  seen  in  Figure  7b.  The  initial 
decay  is  slower  than  for  an  exponential.  This  initial  decay  region  is  also 
usually  quite  ragged  for  many  possible  reasons.  The  interference  effects  of  the 
gage  itself  alters  the  shock  rise.  There  are  reflections  off  the  gage  mount, 
tabs,  and  feedthrough.  Probably  the  most  common  reason  for  departure  at  the 
peak  is  that  the  charges  tested  are  rarely  centrally  initiated  spheres.  They 
are  usually  cylinders  with  various  length  to  diameter  ratios;  some  are  cased 
charges.  Edge  effects,  perturbations,  and  lack  of  symmetry  are  to  be  expected. 
These  signatures  persist  for  distances  far  greater  than  the  normally  believed 
"several  charge  radii". Much  of  what  the  gage  actually  "reports"  is  what  it 
truly  "sees." 

5.2  RECORD  "CORRECTION"  TECHNIQUES 

A  commonly  used,  and  historical  method  to  "correct"  the  shock  wave  for  loss 
in  peak  due  to  the  gage  response  was  to  extrapolate  from  about  one  0  back  to 
one-half  the  transit  time  of  the  gage,  the  zero  reference  for  time  and  distance 
measurements.  When  this  method  was  first  instituted  gages  tended  to  be  larger, 
were  generally  wax  coated,  and  the  gage  rises  were  slower,  rounded  over,  and 
smoothed;  the  gage  filtered  high  frequency  perturbations.  Today,  with  smaller 
gages  with  more  "transparent"  coatings  or  coverings  and  higher  frequency 


^Cole,  Underwater  Explosions. 
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response  instrumentation,  this  extrapolation  technique  becomes  a  uniform  method 
for  idealizing  the  real  or  actual  shock  wave  into  an  exponential.  The  method  is 
still  valid  and  may  well  be  the  only  reasonable  way  to  handle  these  waveforms 
for  comparison  purposes.  This  method,  however,  is  more  a  uniform  idealization 
that  includes  correction  for  the  gage  when  necessary.  Often,  the  gage  size  is 
small  enough  that  the  gage  response  should  not  alter  the  peak  pressure 
measurement.  In  these  cases  the  actual  extrapolated  peak  pressures  are  often 
found  to  be  somewhat  different  than  the  gage  accuracy  predictive  techniques  of 
Chapters  3  and  4.  The  differences  can  be  higher  or  lower  than  the  apparent 
peaks  and  are  generally  within  10  percent  of  the  predicted.  This  departure  from 
the  predicted  is  dependent  upon  the  deviation  of  the  "real"  shock  wave  from  an 
exponential  as  discussed  in  Section  5.1.  No  technique  should  be  used  beyond 
about  a  15  percent  correction. 

This  extrapolation  method  of  idealizing  the  shock  wave  for  uniform  computer 
processing  is  expedient  and  works  well  for  many  types  of  explosive  comparison 
tests;  it  can  mask  detail  in  gage,  model,  or  special  charge  configuration  tests. 
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CHAPTER  6 
CONCLUSIONS 


The  oil-booted  tourmaline  gage  exhibits  most  all  the  ideal  qualities  for  a 
gage.  The  oil-boot  covering  does  not  change  the  rise  response  shape  or  transfer 
function;  the  gage  just  appears  to  be  larger  (about  1.5  times). 

The  technique  of  expanded  rise  comparisons  between  the  real  and  geometric 
is  a  valid  and  useful  method  to  evaluate  gage  coating  materials.  This  along 
with  gage  calibration  constant  (KA)  comparisons  describes  the  characteristics  of 
the  real  gage. 

Accuracy  predictions  for  selecting  a  gage  size  are  made  by  using  the 
predicted  ratio  8/tj)  along  with  Figures  3  and  4  and/or  Table  1.  The  actual 
transit  time  can  be  used  if  known  or  a  good  approximation  of  1.5  times  the 
geometric  rise  can  be  used.  The  shock  duration  6  is  pedicted  from  similitude 
equations. 

A  good  uniform  technique  to  idealize  the  recorded  shock  wave  to  an 
exponential  for  analysis  comparison  purposes  is  to  extrapolate  the  waveform  from 
about  one  6  back  to  1/2  the  transit  time  tj).  Should  a  gage  correction  be 
necessary,  this  technique  corrects  for  it  as  well. 

The  geometric  predicted  peak  may  vary  from  the  apparant  or  the  extrapolated 
peak  for  the  many  reasons  given  in  Chapter  5.  The  differences  are  primarily 
dependent  upon  the  extent  to  which  the  shock  wave  approaches  a  true 
exponential.  These  differences  should  be  within  10  percent.  No  technique 
should  be  considered  valid  if  the  differences  exceed  15  percent. 
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APPENDIX  A 

GAGE  RESPONSE  DERIVATION 


The  tourmaline  gage  response  is  derived  by  treating  the  gage  as  a  circular 
disc  facing  an  incoming  shock  wave  edge-on.  The  disc  with  its  labeled  axes 
appears  below  in  Figure  A-l: 


SHOCK  SPEED  -  c 
2R/c  * 


FIGURE  A-1  CIRCULAR  GAGE  ELEMENT 

The  equation  for  the  circle  shown  in  Figure  A-l  is: 

(t  -  R/c)2  +  p2  =  (R/c)2 

P  =  ((IM  -  t)2)1/2 
c 

The  step  response  g(t)  is  a  function  of  the  area  under  the  curve: 

rt  t 

g(t)  -  2)  P  dt  *  2j  ((25JL  -  t)2)1/2  dt  (3) 

0  0  c 

It  is  convenient  to  normalize  this  to  get  the  unit  step  response.  This  is  done 
by  dividing  by  the  total  area  of  the  circle: 

g(t)N  =  - 2 - f  P  dt  (4) 

ir(R/c)2*u 

The  impulse  response  is  simply  the  differential  of  the  step  response,  thus: 

g'(t)  -  i.  Cg(t)K)  -  - 2 -  (i^l-t2)1/2  (5) 

w(R/c)2  c 

A-L 


(1) 

(2) 
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The  response  of  the  circular  gage  element  or  disc  to  an  exponential  shock 
wave  is  obtained  by  convolving  the  impulse  response  with  the  exponential  of  the 
form: 


P  =  Pm  e_t/e  (6) 

where  8  is  the  shock  wave  duration  or  time  constant,  or  point  where  P  has  the 
value  Pm/e  (e  =  2.718).  The  peak  Pm  is  normalized  to  one  to  get  simply: 

P  =  1  •  e-t/0  (7) 

The  form  of  the  convolution  integral  appears  below: 

C(t)  =/  g'(X)  P(t  -  X)  dX  (8) 

0 

The  expression  for  the  input,  the  exponential  shock  wave,  must  be  convolved  with 
the  unit  impulse  response  to  get  the  response  of  the  circular  gage  element: 

C(t)  =  2  f (2R  X  -  x2)l/2  e"(t  ~  x)/6  dX  (9) 

n(R/c)2  J  c 

A  Fortran  program  was  written  and  the  above  equation  was  evaluated  for  the  rise 
of  the  normalized  gage  response  for  many  values  of  8,  with  a  fixed  value  for 
tj),  the  transit  time  of  the  gage  or  the  time  to  cross  the  diameter  of  the 
element: 


tD  *  2R/c  (10) 

Fifty  points  across  the  circular  element  were  evaluated  for  each  value  of  8. 

The  transit  time  for  a  1/4-inch  gage  was  used  for  computational  purposes.  Final 
results  were  tabulated  for  the  ratio  of  ©/ tjj  to  reduce  and  normalize  these 
so  that  they  would  be  valid  without  regard  to  gage  size. 

It  should  be  noted  that  all  amplitudes  were  zero  before  zero  time,  the  time 
at  which  the  shock  wave  reaches  the  gage  edge.  The  unit  step  and  the 
exponential  jump  instantaneously  from  0  to  their  value  of  1  at  zero  time. 
Integrations  were  taken  across  the  gage  edge-to-edge ,  or  from  zero  to  2R/c.  The 
impulse  response  has  value  zero  before  zero  time,  and  after  the  transit  time  of 
the  gage  2R/c.  Results  were  both  tabulated  and  plotted  as  in  Figure  2. 


A- 2 
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